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Connecting the Kontsevich-Witten and Hodge 
tau-functions by the GL( oo) operators 
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Abstract 

In this paper, we present an explicit formula that connects the Kontsevich- 
Witten tau-function and the Hodge tau-function by differential operators belong¬ 
ing to the GL{ oo) group. Indeed, we show that the two tau-functions can be 
connected using Virasoro operators. This proves a conjecture posted by Alexan¬ 
drov in [Ij. 

1 Introduction 

The tau-functions have been playing a very important role in the study of inte- 
grable systems. One of the known facts is that the set of all tau-functions of the 
KP hierarchy forms an orbit under the action of the group GL{ oo). This group 

is associated to the infinite dimensional Lie algebra gl(oo), which is a central ex¬ 
tension of the fll(oo) algebra. Many well-known tau-functions arose in algebraic 
geometry, such as the Kontsevich-Witten tau-function, Hodge tau-function and 
Hurwitz tau-function. Since the group GL{ oo) preserves the KP integrability, re¬ 
lations between such tau-functions using operators of this type are particularly 
interesting. For example, Kazarian proved in [14] that the Hodge tau-function 
exp (Fj{(u,q)) is indeed a tau-function for the KP hierarchy, where Fn(u,q) is 
obtained from the Hodge series Fu{u,t) after a change of variables derived in 
[Tlj . The main idea of his proof was to show that expg)) is connected 
to the Hurwitz tau-function exp (H) by a GL{ oo) operator, where H is the gen¬ 
erating function of simple Hurwitz numbers and it was known that exp(iL) is a 
tau-function of the KP hierarchy (cf. |23j.[15j). 
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The Kontsevich-Witten tau-function, denoted by exp(Fx ((/)), satisfies the 
KdV hierarchy (cf. [16]). The function Fx(q) is obtained from the generating 
function Fx(t) for the intersection numbers of ?/>-classes after the variable change 
ti = (2i— l)\\q 2 i+\. Using Mumford’s theorem in [22], one can transform Fx(t) into 
the Hodge series Fff{u,t) using a differential operator W described by equation 
m, namely, 

exp (F H (u,t)) = e 117 • exp (F K (t)). (1) 

However, W does not belong to the gl(oo) algebra. Recently, Alexandrov posted 
a conjecture in [I] stating that there exists a G'L(oo) operator connecting the 
Kontsevich-Witten tau-function exp (Fxiq)) and Hurwitz tau-function exp (H). 
This operator consists of two parts, where the first part connecting exp (H) and 
exp (Fi{(u,q)) is already know in [13]. Note that the tau-function exp (H) has a 
description of the cut-and-join operator el which belongs to the GL{ oo) group. 
Hence, if this conjecture is proved, not only will it allow us to have a GL{ oo) op¬ 
erator description of the Kontsevich-Witten tau-function, but also provides a pos¬ 
sible way to derive the Virasoro constraints for the Hurwitz tau-function. Later, 
Alexandrov refined this conjecture as Conjecture 2.1 in [2] stating that there exists 
a GL{ oo) operator G + formed by Virasoro operators, such that 

G + • exp (F H (u,q)) = exp (F K (q)). 


He also verified this conjecture up-to an unknown factor, which is a Taylor series 
in variable u, and hence confirmed a weaker version of this conjecture (see Section 

4)- 

In this paper, we give a complete proof of Alexandrov’s conjecture in [Ij. Our 
method is totally different from the method used in |2|. In our proof, we can see 
a clear relation between actions of the differential operator W and the Virasoro 
operators. We will also give a more explicit way for computing coefficients of the 
Virasoro operators. We first present an explicit GL{ oo) operator that connects 
the two tau-functions exp (Fji(u,q)) and exp (Fx(q))'- 

Theorem 1 The relation between the two tau-functions exp (Fh(u, q)) and exp (Fx(q)) 
can be written as the following: 


exp (F h (u, q)) = exp (^ a m u' n L m ) exp(P) • exp (. F K (q )), 


m> 0 

where L m ,m > 1, is the Virasoro operator 

L m = ^ (k + m) qt - si — + - ^ ab^- 

k>0,k-\-m>0 a-\-b=m 

The coefficients {a m } can be computed from the equation 
exp(^ • 2 = (—2 log(l - 


m> 0 


dz' 


l+Z' 1 + Z 


( 2 ) 


(3) 
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And P is a first order differential operator of the form 


P = ~Y^ b 2k+lU 2k 

k =1 


d 

dq2k+3 ’ 


where the numbers {62/c+i} are uniquely determined by the recursion relation, 

n— 1 

(n + 1 )b n = 1 - ^ kb k b n+ i_ fc (4) 

k=2 

with b\ = 1 ,62 = 1/3. 

The coefficients {a m } are first determined by Lemma[9]below. In Section 3, we 
will prove this main result. The method is to decompose the two operators exp(IT) 
and exp (Yl m >0 a mF m L m ) into several factors using the Zassenhaus formula, then 
study the action of these factors and compare them. One of the difficulties in the 
proof lies on the fact that Fx(q) only involves odd variables < 722+1 while Fji(u,q) 
involves all variables qi. This causes the problem that the decompositions of 
the above two operators do not exactly match with each other. We will show 
that the differences of these operators do not matter when acting on the function 
exp (F^iq))- For this purpose, we need to implement an integral transform on 
functions using Gaussian integrals. We also need to use properties of the gamma 
and beta functions. More details can be found in Section 3.3. 

We can also use an explicit function h(z) defined in Corollary 1121 to determine 
the coefficients {a m }. The power series h and the function on the right hand side 
of equation Q are both related to the Lambert W function, about which we will 
explain more in Section 2.3. After a brief introduction in Section 2.2, we can 
see that the Virasoro operators and P all belong to the infinite dimensional Lie 
algebra gl(oo) acting on the space of KP solutions. Therefore the KP integrability 
is preserved in our expression. Note that the operator exp(P) performs a shift 
of arguments on a series. The extra parameter u in P guarantees that exp(P) • 
exp (Fxiq)) is a well-defined formal series. Furthermore, due to the fact that the 
Kontsevich-Witten tau-function satisfies the Virasoro constraints, an application 
of Theorem [T] is the following 

Corollary 2 There exists a sequence of numbers {e m }, such that 


exp {F h (u, q)) = exp e m u m L m ) • exp (F K (q)). 

m>0 

The coefficients {e m } are determined by equation f W\ ) below. 

This gives a complete proof of Alexandrov’s conjecture in pLJ . We have computed 
the first few terms of {e m } by hand, and they coincide with Alexandrov’s data. We 
will explain this in Section 4. Moreover, since the Kontsevich-Witten tau-function 
satisfies the KdV hierarchy, Theorem [T] gives an alternative proof for Kazarian’s 
theorem that exp {Fh(u, q)) is a tau-function of the KP hierarchy. This proof does 
not need the theory of Hurwitz numbers and ELSV formula. 
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2 Preliminaries 


2.1 Kontsevich-Witten and Hodge series 

Let M g ^ n be the moduli space of complex stable curves of genus g with n marked 
points, and ipi be the first Chern class of the cotangent space over M g n at the ith 
marked point. The intersections of the ^-classes are evaluated by the integral: 

< r dl ...T dn >= [_ 

”Mg } n 

The Kontsevich-Witten generating function in variables t k is defined as 


F r<(t) = Y] < 


k 0 k\ 

W 


> 


j .ko jfci 
L 0 L 1 

ko'.h'. 


Let A j be the jth Chern class of the Hodge bundle over M g ^ n whose fibers over each 
curve is the space of holomorphic one-forms on that curve. The Hodge integrals 
are the intersection numbers of the form 

< A jT dl ...T dn >= [_ A rff 1 ...i l>i n . 

JMg, n 

They are defined to be zero when the numbers j and di do not satisfy the condition 


n 

j + ^ dj = dim(M gin ) = 2>g - 3 + n. 

i=l 


Hodge series is defined as 

F H (u,t) = J](-l ) J < A• • • > • • • 

where u is the parameter marking the A-class. In fact, Fff(0,t) = Fx(t). As 
we have mentioned before, the two functions Fx(t) and Fn(u,t) are related by 
operator exp (IK) as in equation (P), (cf.[8].[T0]). where 




d 


k> 1 


2k{2k - 1) y dt 2k 


j>0 


d 1 
dti + 2k-l 2 


E (- 1 )' 


d 2 


i-\-j=2k—2 


dtidtj ’ 


(5) 
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Here B 2 k is the Bernoulli numbers defined by: 


j. 00 j m 

z — J ml 


e L — 1 


ra =0 


Note that VP does not belong to the gl(oo) algebra. 

Now let 

FK (?) = i? *'( t )lt fc = ( 2 fe-l)!!q 2fc+1 ' 

It is well-known that exp (Fj{(q)) is a tau-function for the KdV hierarchy jl~ 6 j . 
And exp (Fff(u,t)) is a tau-function for the KP hierarchy after the a change of 
variables derived in |T3]. It can be described in the following way. Let 

D = (u + z) 2 z ( 6 ) 
Consider the following sequence of polynomials cpk(u,z ): 


2fc+l 

<j>o(u,z) = z , 4>k{u,z) = D k z = u 2k+l ~ j z j , 

3 = 1 

for some constants , where a^t+i = — 1)!!- Let 4>k(u,q ) be the polynomial 

in variables % and u obtained by replacing z m with q m in 4>k(u, z). Then the 
variable change is in the form 

2fc+l 

tk = <t>k{u,q) = ^2 u 2k+1 ~ 3 qj, 

3 = 1 

For example, 

to = qi 

t\ = u 2 qi + 2 uq 2 + q 3 

t 2 = u A q\ + 6 u 3 q 2 + 12 u 2 q 3 + 10u ^ 4 + 3 q 5 

£3 = u 6 qi + 14u 5 (/2 + 61u 4 (?3 + 124it 3 (/4 + lSlu 2 ^ + 70uqe + lhq^ 


Note that, by definition of the double factorial, 0!! = (—1)!! = 1. Let 

F H (u,q) = F H (u,t)\ tk= ^ q y 

If we set u = 0, the polynomial ct>k(u,q ) only has the term (2k — l)!!<Z 2 fc+ij and 
exp(Ftf(0, q)) becomes the Kontsevich-Witten tau-function exp (Fx(q))- 
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2.2 The transformation group GL{ oo) 

The space of tau-functions for the KP hierarchy forms a group orbit [7j. We denote 
this group by GL( oo), and it is defined from the infinite dimensional Lie algebra 
gl(oo) via the exponential map. Note that, strictly speaking, the “group elements” 
in GL{ oo) constructed in this way do not belong to a well-defined group, because 
the products of these elements might be divergent. Hence when we want to use 
a “group element” formed by a product of other such elements, we need to check 
whether this product is well-defined. 

Here, we introduce some sample operators that belong to the Lie algebra gl(oo) 
and will be used in our results. We refer the readers to [2Cjj or other related articles 
for more details. Suppose the tau-functions are in variables qk . The operators qk 
(i.e. multiplication by qk) and d/dqk belong to gl(oo). It follows that adding a 
linear function to the solution or a shift of arguments preserve the solutions of the 
KP hierarchy. The Virasoro operators {L m } that generate the Virasoro algebra 
with central charge c = 1 also belong to this Lie algebra. For m € Z, setting 
qj = 0 and d/dqj = 0 for j < 1, we can write L m as 


Lm — 


\ ^ . . d 1 \ ^ cP' 1 \ ^ 


k>0,k-\-m>0 


a+b=m 


dq a dqb 


qaQb, 


a-\-b=—m 


and they satisfy the commutator relation 

] 

[L m , L n ] — ( Tfl 7l)L m _j_ n -j- ——Sm^-n^TTl Tfl). 

In particular, the degree operator Lq = kqkd/dqk■ In Section 3, we will discuss 
how we deal with the exponential of Virasoro operators when they act on functions. 


2.3 Gamma function and beta function 

The gamma function T(z) is an extension of the factorial function. If z is a 
positive integer n, then 

T(n + 1) = n\. 

And when z is a complex number with positive real part, then the following 
equation holds 

T(z + 1) = zT{z). 

It can be defined by a definite integral, known as the Euler integral of the 
second kind, 

COO 

r(z)= / x^e^dx. (7) 

Jo 

This definition is only valid for complex number 2 with positive real part. The 
function r(T) can be extended to a meromorphic function on the whole complex 
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plane with simple poles at non-positive integers. Furthermore, using Stirling’s 
approximation, one can obtain an asymptotic expansion for T(z) with the value 
of \z\ large enough: 

OO 

r(z) ~ z z -h- z V 2 ^^CiZ-\ (8) 

i =0 

with Co = 1. The coefficients {C*} play an important role here in this paper. In 
fact, 

Ci = (2i + 1)!!6 2 *+i, (9) 

where the coefficients { 624 + 1 } appear in the relation dH) (cf. [[18]). If we take the 
logarithm on r(z), the coefficients B 2 k appearing in the expansion are simply the 
Bernoulli numbers: 


where 


log(r(z)) ~ (z - 5 ) log(z) - z + \ log( 2 vr) + <B(z), 


»(*) = E 


-^2 k 2 fc+l 


2k(2k - 1 )' 


This implies 

OO 

e®( z ) = E CiZ^ 1 . (10) 

4=0 

If we expand the left hand side of the above equation, we obtain the general 
formula as 


= 1 + 


E 

m> 0 


ml 


E 


n 


B 


2k n 


k\ ,...km> 1 j —1 


2kj{2kj - 1 ) 


y m ETLi 


Hence 


Ci = E 

m= 1 


1 

m\ 


e n 

Ej=i 2 fcj=i+mj=l 

All 1 


B Zkj 

2kj (2kj - 1 )' 


( 11 ) 


The fact that %$(z) is an odd function immediately implies the following relation 
on the coefficients {C,},i > 0. 

Lemma 3 For all k > 1, we have 


k 

^(- 1 )^ = 0 . 

4=0 

Furthermore, for i > 1, C,; has some nice combinatorial interpretations. It can be 
explicitly given by 

r _ y^/ ,sk dz{2i + 2k,k) 

* ^ [ ’ 2 *+k(i + k)\ ’ 
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where d^{n,k) is the number of permutations on n elements consisting of k per¬ 
mutation cycles with length greater than 2, (cf. 0). 

The following results from [18] gives us an equation that packs all the coeffi¬ 
cients {bi} into its solution. We first observe that the function ve l ~ v is increasing 
for v € [0,1], and decreasing for v € [1, oo). Now, consider the equation 

ve x - v = e-\ x \ ( 12 ) 


The above equation defines a one-to-one relation between v(x) in [l,oo) and x in 
[0,oo). In this case, v increases with x, and we can write v as a power series of x 
in the form 

OO 

v = 1 T ^2 (13) 

2=1 

The existence of the power series v(x) is guaranteed by the Lagrange inversion 
theorem [T2]- Then, differentiate the equation by x on both sides, we have 

v (v — 1) = xv. (14) 

This gives us exactly the relation ([U) on the coefficients {bi} with b\ = 1. Since 
v > 0, we have b\ = 1. The choice of —1 for b± leads to another solution w(x) for 
equation (fT21) with w € (0, l],x € [0, oo), and w is decreasing with respect to x. 
In this paper, we always set b\ = 1. Then w{x) is written as 

OO 

w = 1 + ^(— l) l biX l . (15) 

2=1 


Note that both v and w satisfy the equation (HD- Also, it was pointed out in 
@1 that the series v (and w) converges for all finite x using the properties of the 
Lambert W function. In fact, many power series appeared in our context are 
related to v and w. 

The Euler integral of the first kind defines the Beta function B(x,y ) for 
complex numbers x and y with positive real parts: 


B(x,y)= f t x -\l-ty- l dt. 
Jo 

And it has many other forms including 


B(x,y) 


£ (x)T(y) 

T(x + y) ' 


(16) 


(17) 


All the formulas and results mentioned in this section will be used in our later 
proofs. 



3 From Fx(q) to Fjj{u , q) 


Before proving Theorem |Tj in this section, we first give an outline for the main 
ideas of the proof. What we aim to show is that after switched to the variables 
qi, equation ([TJ implies equation (f2j), i.e. the following diagram commutes: 


Equation Q 

exp(F k (t )) --► exp(F H (u, t)) 


t k = (2 k — l)!!<?2fc+i 




exp(F k (q)) 


Equation m 



exp(F H (u,q)) 


For this purpose, we want to study the relation between the operators in equation 
© and exp(W) after the change of variables. The crucial steps in our proof can 
be summarized into several propositions, which will be proved later. 

Exponential of second order differential operators are very complicated in gen¬ 
eral. Therefore it is rather difficult to compare operators exp(^ m>0 a m L m ) and 
exp(W) directly. On the other hand, the exponential of a derivation is much sim¬ 
pler since when acting on a function, it just behaves like a change of variables. 
Inspired by this fact, we will decompose operators W and L m into several factors, 
and first compare the exponential of the derivation part of these operators. Such 
comparison will determine the coefficients a m in Theorem [lj 
For the operator W defined in equation ©. Let 

W = !B t + \Qf + P 0 , (18) 

where 


- 2 ^ 01,(01,-U l^ tj 


d 


^ 2k(2k — 1) ^ l dti+2k-\ 


P ° = -E 


B 2k u 2 ^ d 


^ 2k(2k - 1) dt 2k ’ 

pW _ y' B 2k U 2 ( 2k ~ 1 '> y+1 d 2 

^0 - 2 ^ 0b(0.b-1\ 


2k(2k - 1) 

k> 1 v ’ ij> 0 

i-\-j=2k —2 


dtidtj 


(19) 


Using Zassenhaus formula (cf. equation (1371) 1. we can decompose the operator 
exp(W) into the following form : 


Proposition 4 

exp(TU) = exp( s B t )exp(^Q t H ")exp(P f ), 


where Pt is a first order differential operator defined in Lemma 
second order differential operator defined in Lemma 1731 




and QY is a 
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This proposition follows from Lernma ll4l and Lemma [T5l in Section 3.2. Using this 
proposition, we can transform equation © into 

exp (F h (u, t)) = exp(55 t ) exp (^QY) exp (P t ) • exp (F K (t)). (20) 

The operator exp(P) in equation © is obtained immediately from exp(i-)) in 
Lemma [HI after the change of variables t k = [2k — l)!!? 2 fc+i • 

For m > 0, let L m = X m + Y m , where 


Xqm. — 


^2 (k + rn)q k 


d 


k>0,k-\-m>0 




V — 

1 m. — 


E 

a-\-b=m 


ab 


& 2 


dq a dq b 


( 21 ) 


We then compare operators exp( < B i ) and exp(]U m>0 o m X m ), where the coefficients 
a m are chosen in the way such that the following proposition holds: 

Proposition 5 For any polynomial or power series G in variables u and t k , k > 
0, we have 


exp <E a m ^ -X m ) 


m> 0 


t k =(2k-l)\\q 2k+1 


= e 


- U&t 


G 


t k =rp k (u,q) 


where the numbers { a m } are determined by equation ©, and t k = 4>k[u,q ) is 
defined in Section 2.1 before. 

We need to explain the notation “ • ” used in the above formula. For a differential 
operator T> and a function /, we always use the notation “ T> ■ f ” to represent 
the action of V on /, in order to distinguish from “ T>f ” which could represent 
the multiplication of T> and / as operators when / is considered as an operator 
acting on the space of functions by multiplication. 

We will prove the above proposition using Lemma 0 which is how we obtain 
the series on the right hand side of © to represent the numbers {a m } in the first 
place. 

Let 

n w = Q W I 

^9 b/i lt fc =(2fe-l)!! 92fc+1 ' 

Applying the change of variables t k = </>fc(u, q) on both sides of equation (f20l) . and 
using Proposition [5] for G = exp {^Q] v ) exp (P t ) ■ exp(Fx(^)), we can establish the 
following connection between exp(F#(u, q)) and exp(F/i ((/)): 

Corollary 6 

exp (F H [u,q)) = exp( E a m u m X m ) expexp (P) ■ exp (F K (q)). 

m> 0 
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This expression is very close to equation (J3J) in Theorem [T] now. To complete the 
proof of Theorem [TJ, we need to compare the following two operators: 

exp( ^ a m u m X m ) exp (^Q^) and exp( ^ a m u m L m ). 

m> 0 m> 0 

Since L m = X m + Y m , we can use Zassenhaus formula to obtain the following 
decomposition 

Proposition 7 

exp(^ a m u m L m ) = exp(^ a rn u m X rn ) exp(-Q + ), 

m >0 m >0 

where Q + is a second order differential operator defined in Corollary HT1 

It turns out that the operators Q + and Q ^ are not the same. This fact explains in 
some sense why it is difficult to compare operators exp(^ m>0 a m L m ) exp(P) and 
exp(W) directly without decomposing them first. The application of Zassenhaus 
formula in our proof singles out the major difference between these operators. 
Fortunately such difference does not affect Theorem [T] due to the following propo¬ 
sition: 

Proposition 8 Let Q+ dd he the sum of terms in Q + involving only odd variables 
q 2 k+i- Then, 

O w — o + 

We would like to remark that the expressions for Q + and Q^ which arose during 
the process of applying Zassenhaus formula are very different. The relation of 
these two operators as given in the above proposition is not obvious at all when 
we first saw them. To prove this proposition, we need to express these operators in 
simpler forms. Hence in Section 3.2, we introduce the formal power series Q B (x , y ) 
in Lemma fT5l and Q(x,y) in equation (ITTIi to represent QY and Q + respectively. 
We then perform a transformation of Q + using the Gaussian integral and also 
use properties of Gamma function and beta functions to prove Proposition [8] in 
Section 3.3. 

Combining the above results, we have 
exp (F H (u,q)) 

= ex P (^2 a ™u m X m ) exp(^Q+ d ) exp(P) • exp ( F K (q )) 

m> 0 

= exp ( ^2 a m u rn X m ) exp(^Q + ) exp(P) • exp ( F K (q )) 
m>0 

= exp ( ^2 a m u m L m ) exp(P) • exp ( F K (q )). 

m> 0 
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The first equality above follows from Corollary [6] and Proposition [HJ The sec¬ 
ond equlity follows from the fact that exp (Fx(q)) and the operator exp(P) only 
contain odd variables q^fc+i- The third equality follows from Proposition 0 This 
completes the proof of Theorem |TJ As mentioned before, we can replace the oper¬ 
ator exp(P) with the exponential of a linear combinations of Virasoro operators, 
using the fact that the Kontsevich-Witten tau-function satisfies the Virasoro con¬ 
straints. This leads to a proof of Corollary 2, which will be given in Section 
4. 

In the rest part of this section, we will prove the above propositions. 


3.1 Proof Proposition [5] 


In this subsection, will prove Proposition [5] and also discuss properties of the 
sequence of numbers a m appeared in this proposition. We first discuss change of 
variables using differential operators in the form e x , where X is a derivation. This 
is due to the following basic fact. Suppose we have two elements / and g from a 
commutative algebra, and X is a derivation on the algebra satisfying Leibniz rule, 
that is, X ■ ( fg ) = (X ■ f)g + f(X ■ g). Then, it is easy to verify, by induction for 
instance, that 


x" •(/«)= E (i+i)!—r 


X'-l xi -g 


i-\-j=n 


and 

e X ■ (fg) = (e x ■ f)(e x ■ g). 


It follows that applying an exponential of a derivation on a polynomial function or 
a series simply performs a change of variables. This works for series with infinite 
number of variables too. 

In this paper, many problems involve differential operators and series with 
infinitely many variables. However, we will reduce such problems to those with 
just one variable and containing infinitely many terms. This method can simplify 
the computation process, and it has been frequently used in many papers, e.g. [131 . 
In other words, we will use power series in one variable to describe the action of 
the operators used in our context. For the knowledge of power series, including the 
compositional and multiplicative inverse of a power series, derivatives, composition 
of two power series and products, we refer the readers to the book }12j . 

First, we will consider X m and T>t, which are derivations appeared in the Vi¬ 
rasoro operator L m and the operator W, and are defined in equations (12T1) and (HU 
respectively. Since we need to deal with operators in the form exp(^ m>0 a m u rn L rn ) 
and exp(IF), it is necessary to discuss how the two operators 


exp( E a m v m X m ) and exp (23i) 

m> 0 
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behave as changes of variables. The connection between these operators is given 
by Proposition [5j Since X m and are derivations, to prove Proposition [5l we 
only need to consider the case G(t) = t n , which has the following form 

Lemma 9 There exists a unique sequence of numbers > 1, such that, 

for all n > 0, 


exp <E a m U • <72n+l 

m> 0 


1 

(2tt. — 1)!! 





To prove this lemma, we will use operators with one variable z to help us with 
the calculation. First we introduce an isomorphism ©i between the vector space 
of formal power series in variable z and the vector space of linear functions in 
variables qf. 


0 i : (E 


a{Z l | are constants} —> anqi | a* are constants} 

i> 1 i> 1 


@i C^atiZ 1 ) = y^ y aiqi. 


i> 1 


i> 1 


For a sequence of numbers {a k }, let 


= E a ^ 1+k 

k =1 


d_ 

dz ’ 


*J = E 


= 2^ a k z 
k =1 


i-k_d 
dz' 


The action of e® z on z n for n E Z is of the form: 


m\ 


m> 0 


f r - 2 ” 


i 


m—1 


=z™ + ^ ^ na fcl (n ± fci)a fc2 ... (n ± ^ 

m>0 fci,...fc m >l j=l 


^±(Er=i fc p. 


( 22 ) 


The result of this action is a formal Laurent series. We can use this series to obtain 
the action of the following two operators on variables q n for n > 1 respectively: 

ex P(X/ UmX rn ), exp( E 

m<0 m>0 

For example, for the second case above, we have the following formula: 


exp( E ‘ Qn 

m> 0 


m—1 


=9n + 2^ 2^ nafc i ( n “ kl ) flfc 2 • • • ( n - XI 9n-(E7=i fcj) j 

m>0 ki,...k m >l j=1 
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and this gives us a linear function in variables (fr. The coefficient of qi on the right 
hand side of the above equation is the same as the coefficient of z l with i > 1, in 
the series exp(<f>7) • z n . More precisely, let 

OO 

F{z) = e** -z n = z n + Y^ A\ n) z n -\ (23) 

i =1 


and define ( F(z)) + to be 


n— 1 

(F(2)) + =2» + £4 < "V , -‘- 

l=l 

Then, 

0i ((e #2 • z n )+) = exp(^ a m X m ) ■ q n . 

m> 0 

Note that in this formula, we can not replace z n by a formal power series p(z) = 
ot n z n since both sides of this equation are not well-defined in this case. 
The action of exp(^) m>0 a m u rn X m ) on q n can be obtained by the variable change 
qk —> u~ k qk in the above equation, that is 


exp( E dm u m X n 


m> 0 

Similarly, we can deduce that 


n— 1 

q n = q n + 

i= 1 


+ ^2 A^u'qn-i. 


0! • z n j = exp(^ a- m X m ) ■ q n . 


m< 0 


More generally, for any power series p{z) = X^n>i a nZ n , we have 
01 ( -p{z2j = exp ('^2 a- m X m ) • 0i(p(z)). 


m< 0 


On the other hand, the action of exp(iBf) on t n is of the form: 
exp(«B t ) • t n 

=‘»+e^ e n 


rn >0 *+m=E7Li 2k j j 


1-L 2kj{2kj - 1) 


(24) 


(25) 




i =0 
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where the last step is implied by equation (fTTT) . Applying the variable change 


2k+l 

tk = 4>k(u,q) = 22 a f )u2k+1 ~ J qj, 


3 = 1 


to the last equation, we have 



tk=<Pk{u,q) 


n 

= 22 C iU 2l (l) n -i(u,q) 

i =0 

n 2(n—i)+l 

=E C < E 

*=o j =i 


Hence, in order to prove Lemma 0 we only need to consider the case when u = 1, 
that is, 

1 n — 

ex P( a mXm ) ' 92n+l = 7^^ :jTjj ^ ^ Cj(j) n —i{\, q). (26) 

m>0 ' ' i=0 

The general case can also be obtained from the above equation after the change 
q k ->• u~ k q k . Let 

D = {1 + z)2z lL ’ 

which is the operator D defined by equation © with u = 1. Under the isomor¬ 
phism @1 we introduced before, the left hand side of equation (1261) is equal to 
0i • z 2n+1 )j2j • And on the right hand side, by the definition of 4> k (u, z) and 

4> k (u,q) introduced in Section 2.1, we can see that 

0i ( D n ~ l ■ z) = ©i (0 n _j(l,z)) = (j) n -i(l,q). 

Therefore Lemma O follows from the next lemma: 

Lemma 10 There exists a unique sequence of numbers {a m },m > 1, for the 
operator such that, for all n > 0, 

1 n 

(e * r ■ " 2 ” +1)+ = (27mTT E c- 0 ” - '' P?) 

v i =0 

Proof: We hrst prove the uniqueness part by induction. Assume there exists a 
sequence {a m },m > 1, for operator <3?“, such that equation (1271) holds for all 
n > 0. When n = 0, both sides of equation (1271) are equal to For n = 1, we 
have, on the left hand side of equation (1271) . 

{e^ z ■ z 3 ) + = z 3 + 3aiz 2 + (3 af + 3ci2 )z. 
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On the right hand side, we have 

13 

Dz + C\ z = z 3 + 2 z 2 + —. 

This gives us the only solution a\ = 2/3, a 2 = —1/12. Now we assume {dj}, 
for 1 < i < k, are all uniquely determined by equation (1271) . Then, we choose 
n large enough such that 2 n > k. By equation (1221) . the coefficient of z 2n ~ k in 
(e®* . z 2 n+l )+ ig 


^ — (2n + l)afc+i + 

^ m—1 

X] ( 2n + 1 ) a ^i ■ ■ ■ (2n + 1 - X] %)a fem 2 : 2n+1_(E ^i kj) . 

m>2 m ' YJ'h kj=k+l j =1 

The summation in the second term on the right hand side of the above equation 
consists of only {ai,..., a*,}, and 1 is uniquely given by the right hand side 

of equation (1271) . Hence there is only one solution for a^+i- This proves the 
uniqueness. 

Now we discuss the existence part. In order to find the coefficients {a m }, we 
need to find a series / such that / = exp(<hj) • 2 . To do so, we first define a 
function f(z) to be 


f(z) = (—2 log(l - 


1 


1 + z' 


Note that / is a solution of the equation 

Df = f- 


(28) 


(29) 


And its asymptotic expansion, which we also denote by /, can be computed in 
the following way (we explain why choosing this function in the remark later). By 
the equation above, we have 


/ 

2 


2 

2 


~2 Z 


1 2 


log(l + 

4*- 2 

1o &(i+5Ft) + l^pr 

— 1-2 2 

2 Z 


E^=i(-i)^-Er=i(-i)^- n 

i 

°° _ 1 

1 + 2 V(-l) n_1 -- z 2 ~ n + ■■■■ 

2 -—' n 


n =3 
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Hence the expansion of / is of the form 


f=Z+ 3~U Z 


-1 


Note that the coefficient of z in / is set to be 1, not —1. The above series / 
determines a unique set of coefficients {a m } for <I>J. To see this, we set / = 
exp(<3?j) • z, then by equation (1231) . for n = 1 and k > 1, 


/ = 5> S‘V 


where, by equation (1221) . 

= oi; 

k —1 ^ m— 1 

A k )=a k+^2^ a ki(l ~ ... (1 - kjhkm- 

m=2 ’ i k j= k 3=1 

kj> 1 

Using induction on k, we can see that a fixed set {Hj 1 '*} uniquely determines {a 
and wfce versa. From now on, we always assume { a m } are determined by 

ex P(3T ) -z = f, (30) 

where / is the function (|28D . The existence part of Lemma flOl follows from Lemma 
ei below. In other words, the coefficients {a m } determined here are exactly what 
we need for Lemma ITOl 


□ 

Remark: (1) The series exp(<L“) • z _1 can be seen as the multiplicative inverse 
f~ l of /. Since exp(<f>7) • 1 = 1, we have 

exp(<h;) • (. z~ x z) = (exp(d»-) • z~ l )f = 1. 

The standard definition of the multiplicative inverse of a power series usually 
requires that the series has a constant term (cf. |T2]) ■ However, in our case, the 
series j/z can be seen as a power series in z~ 1 with the constant term 1. Hence 
/ -1 can be understood as z~ 1 (f /z)~ l . From the discussion at the beginning of 
this section, we can deduce that, for k € Z, 


and for negative integer k, f k is a formal Laurent series. 
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(2) To see where the function / defined in equation (|28ji comes from, we present 
a heuristic argument here. From the previous remark, we can see that equation 
(1271) can be transformed into 


(/ 


2n+l\ 


n 

— Y.C’ d ' 

i =0 


^ (2n — 1)!! 

This gives us a hint that the series f 2n+l behaves like the formal expression 


(2 n- 1)!! 


E C ‘ D ' 


z. 


i =0 


This expression is not well-defined because it contains negative power of operator 
D. However, it gives us a possible equation that / might satisfy. In particular, 
for n = 2, since 


\Y j C i D 2 ~ i -z = \DY J C i D 1 - i -z, 


i =0 


i=0 


we might expect / to satisfy 


which gives us 


f = -D-f = f(D-f), 


f=D-f. 


Solving this ODE we obtain a solution f(z) which is the function defined by 
equation 


□ 


Lemma 11 For all n > 0, 

( / 2n+ 1 )+ 


1 

(2 n- 1)!! 


n 


CiD n - 1 


■ z. 


Proof: By the definition of function / given in equation (1281) . we have 

-if- 2 1 \ — 

e 2 J = (1-) e i+z. 

1 + z J 

Setting (1 + z) -1 = 1 — v and / = x~ l in the above equation leads us to equation 
(HU). Since equation (fT2l) has two solutions v(x) and w(x ), we obtain 

OO OO 

(i +^)- 1 = ±y b 2 l+l r 2i ~ i - y b 2 i r 2i . 

2=0 2=1 
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Now apply D on both sides of the above equation. Since Df = / 3 , we have 
Df~ l = Hence 

OO OO 

2 ; = ±^T(2i + l)& 2 i+i/~ 2 * +1 - ^(2 i)b 2 if~ 2l+2 . 

i=0 i— 1 

It is easy to see that the highest power of z in the series f k = exp(<3?“) • z k is k 
for k £ Z by equation (1221) . Then f k for k < 0 will not contain terms of 2 with 
positive degree. Since the coefficient of z in / is 1, we can only take the positive 
sign in the above equation. This implies that 

OO OO 

f = z- Y( 2 i + l)b 2i+1 f~ 2i+1 + ]T(2i)6 2i r 2i+2 . (31) 

i=l i —1 

From the equation D ■ f = / 3 , we have D n ■ f = (2 n — 1 )!!/ 2n+1 . Then we can 
transform the equation in the lemma into 

n 

(- D n -f) + = Y J C i D n - i -z. (32) 

i=0 

Furthermore, D n ■ f _2 *^ 1 will contain positive degree of z if and only if n > i. 
Since D n ■ f~ 2n is a constant, D n ■ f~ 2t will not contain terms of z with positive 
degree for all n, i > 1. Now we prove the lemma by induction. 

When n = 0, both sides of equation (1321) are equal to z. When n = 1, we apply 
D on equation (13TT) to obtain 

(D • /)+ = (D ■ z + (—36 3 )L> • / _1 )+ 

= D ■ z + C\z, 

where C\ = 363 by the relation @. So equation (FT?!) holds for n = 1. Assume 
equation (1321) holds for all n < k. Let n = k + 1. We apply D k+1 on equation (13T1) 
to obtain 

fc+i 

(D k+1 • /)+ = D k+1 ■ z + ^(-l)(2z + 1 )b 2i+ i(D k+1 ■ f~ 2i+1 )+ 

i=l 
k +1 

= D k+l ■ z + ^(-l) i+1 (2i + 1 )!!6 2m (Z/ +1 -* • /)+ 

i=l 

k -\-1 k-\-l—i 

= D k+1 -z + Y,(-^ i+1 Ci J2 Ck+i-i-jD-i ■ z 

i=l j=0 

k k+l—j 

= D k+i -z+y^(Y1 (-1 ) i+1 QC' fe+1 _ j _o^ i • ^ 

j =0 i= 1 
k 

= D k+l ■z + Y Ck+i-jDi • z, 

j =0 
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where the last step is implied by LemmaO Hence equation (1321) holds for n = k+ 1. 
This completes the proof of the lemma and also completes the proof of Proposition 

□ 


Corollary 12 For the operator with coefficients {a m } determined by equation 
m, we define the series h(z) to be 

h(z) = e* z ■ z. 

For the series w(z) defined by equation m, we have 

. OO 

and 


i =0 


v 7 i =1 


i -2 


In particular, since w(z) is convergent and non-zero, h(z) also converges for all 

Z / OO. 

Proof: We denote the compositional inverse function of / to be ip(z). Since 
z = e~® z e^ z ■ z = e~® z ■ f(z) = f(e~® z ■ z), 

we have 

ip = e~® z ■ z. 


By equation (|22) 1 . we have 
h(z) =e® z ■ z 


m —1 


=Z 


+ E (-!)(-«,)...(-i-E%>< 


ml 

m> 0 




i=i 


On the other hand, 


■0- 1 =e"*‘ • 2 " 1 


m—1 


=2: 


1 + E^T E (-!)(-«,)...(-l-E %)(-«*») 




ml 

m> 0 k\,...k m >l 


3 = 1 


Compare the above two expressions, we obtain 


1 


h{z)=fi> l (z ) and —-r = ip(z ). 

h\z) 
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Since ip is the compositional inverse function of /, by equation (1311) . we have 


■0 = X ibi 


2 -i 


( 33 ) 


i =1 


Since 


we have 


die/dz = X^(—1 YibiZ 1 1 , 

i =1 

oo 

h{ z ) = zQ^i-iy-Hbiz*- 1 )- 1 = - 


i= 1 


d7n/d2 


By the relation w (w — 1) = zw from equation (1141) . we obtain 


h(z) = 


1 


w{z) 


- 1 , 


(34) 


which gives us the form of h(z) used in this corollary and Theorem [lj 


□ 


Note that the coefficients {o m } can be computed easily using either the series f(z) 
or h(z). Also, since h(z) = 1 /w(z) — 1, by equation (fl2l) . we have 


1 


1 + h 


e 1 +'» = e 2 


-4z 2 —l 


(35) 


Corollary 13 Let 


Then 


d 


v( u , z) = exp(— ^2 a m u m z 1+m —) ■ z. 


m= 1 


rj(u, z) = -\ 2 log(l + uz) - 2 + 


1 + uz 


(36) 


Proof: Solving z in terms of h using equation (l35l) . we obtain 


z — \ 2 log(l + h) — 2 + 


1 + h 


When u = 1, the function 77 ( 1 , 2 ) is the compositional inverse function of h{z). 
Therefore 


?7(M) = a / 2 log (1 + z) - 2 + 


1 + ^ 


Note that the coefficient of 2 in 77 ( 1 , 2 ) is set to be 1. It is easy to see that 

7777(77,2) = 7/(1,772). 


This completes the proof. 

□ 
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3.2 Proof of Proposition |4] and Proposition \ 7 \ 

In this subsection, we will prove Proposition [4] and Proposition [7J The main 
ingredient for the proof of these propositions is the Zassenhaus formula, which can 
be considered as the inverse of the Baker-Campbell-Hausdorff (BCH) formula. 
The BCH formula is the expression of Z as 

Z = log(e x e y ) 

for two elements X and Y from some Lie algebra over a field of characteristic zero. 
A general formula was introduced by Eugene Dynkin in [5], expressing Z as a sum 
of nested commutators involving X and Y. If we add another parameter a to Z 
by letting Z = log(e Q ' V e Qi ), and write Z = X^m>o am Z m , the recursion relation 
in [25] allows us to compute Z m starting from Z\ = X + Y. This relation and 
the general formula can be used in principal to construct the series Z to arbitrary 
degree in terms of commutators. However, once we carry out the computation, 
we will meet its limitation very soon, since many identities involving the nested 
commutators, such as the Jacobi identity, are hidden inside the formula. 

We use the notation “ ad? ” to denote the following nested commutator, 

ad° x Y = Y, ad x Y = [X , Y] , ad j x Y = [A, ad j x l Y}. 

The Zassenhaus formula is known to be the inverse (or dual) of the BCH 
formula. It states that we can expand e°T Y +} ) as 

e a{X+Y) = e aX e aP le a?P 2e o?P 3 ( 37 ) 


where, P\ = Y. and for n > 2, P n is a homogeneous Lie polynomial in X and 
Y of degree n, (cf. [17jb The existence of such formula can be shown by using 
induction and the BCH formula on the product 


g -a n P„ _ _ _ e - a p le - a X e a(X+Y)^ 


(38) 


for n > 2. In fact, the result in (cf. [3J) shows that the term P n can be written as 


n —2 


Pn = 


,0 * 

71 * 7 


(- 1 ) 


77,-1 


ad? Y ad n x l - % Y. 


n i\(n — 1 — *)! 

i=0 v ' 


(39) 


The Zassenhaus formula itself is a formal expression involving a product of infinite 
number of elements. When applying this formula on a specific case, one should 
check that the product of elements is well-defined and does not cause any problem 
of divergence. Later we will see that, in our case, the Zassenhaus formula works 
very well. 

We mainly concern with the differential operators from the following three 
types: 


01 = \ x = a 

I i>0j>l 


d 




i+j 


aij are constants 
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02 


i>0 


A 

dti 


/3i are constants 


l 




03 ^ Z 'y ^ 'ya .6 

a,b>0 


d 2 


dt a dt b 


'jab are constants, 7 ab = 7 6a 


Each 0 j is a Lie algebra, and so is the direct sum </i ® < 72 © 93 ■ In particular, 
elements in 02 ® 03 commute with each other. For A, A G 0 i, Y G 02 , Z € 03 , we 
have 


[A, A ] G 01 ; 


E] ^ 1 a i3@i ^ 02i 

[X,Z] = - £ 


a,6>0,j>l 


Suppose we apply the Zassenhaus formula on exp(A + Y). Then we can write the 
operator as 


e x e Pl 


(i) P (i) p(i) 


P 

P 2 


p y 


where, by equation (|39|) . 


pi 1 ) 

± Tl 


(-i r - 1 

n! 




This is because ad\ l Y € 02 for all n > 1, which means that ad l Y ad\ 1 l Y = 0 
for i > 1. Hence all operators PrP commute with each other, and 


OO 

exp (A' + Y) = exp (A) exp(^ P«). (40) 

n= 1 


Also, observe that ad^ 1 Y is a hrst order differential operator that only contains 
operators d/dt a with a > n. Therefore, in the coefficient of d/dt a is 

always finite for all a > 1. Similarly, we can easily obtain 

OO 

exp (A + Z) = exp (A) exp(J^ P (2) ) (41) 

n=l 


with 


pi 2 ) 
J n 


n! 


ad Y ~ 1 Z, 
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( 2 ) 

where P„ ' is a second order differential operator that only contains operators 

/ o\ 

dt a dtb with o + 6 > n — 1. This guarantees that the operator -Fn ; is 

well-defined. Our argument still works if we use variables q 7 with i > 1 for the 
operators. We can also check the validity of the two formulas (1401) and (1411) by 
applying the BCH formula on the right hand sides. And we refer the readers to 
the appendix for more details. 

Now we apply Zassenhaus formula to some exponential operators used in our 
context. We first prove Proposition [I] which is the combination of Lemma 1+41 and 
Lemma [T5l below. We will apply Zassenhaus formula to the decomposion 

w = + \q^ + P 0 

as given in equation (fT 8 l) . It is easy to see that € 0 i, Pq € 02 , and Qq € 03 . 


Lemma 14 


where 


exp(W) = exp(93t + exp (P t ), 


Pt = ~J_2 C t u 2i 
2—1 


d 

dt i+ 1 


(42) 


with Ci defined by equation m>- 

Proof. First we notice that adffif 1 Pq € 02 for all m > 1, and since Po commute 
with Qq V , we have 

5 = ad ** lp °- 

From the discussion before, we can compute Pt in equation (1421) as 


p, = f; LIL 0( c-i Fn . 


m= 1 


ml 


m— 1 

I a ^ t ^0- 


For, m > 1, the formula for the nested commutator adf^ t 1 Pq is in fact 


ad^Po 


B 


2 kn 


^>’”E E 

i=m ki,...,km>l \j = l J 

E7=i+=i 


. 4 z— 2 m 


5 


dtoi- 


2i— m+1 


(43) 


We can check this using induction. In order to simplify the expression of equations 
later, we define the number Bk to be 


Bi. = 


B 


2 k 


2k{2k - 1 )' 


(44) 
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When m = 1 


ad% t P G = P 0 = ~ E B iU A 


Ai-2 


d 


i= 1 


dt- 


2 i 


Assume the formula (143h holds for m = n,n > 1. Then, when m = n + 1, we have 

adk t P 0 = [®t,ad57 1 P 0 ] ■ 

For k n+ 1 > 1, since 


E p 


3 


9 


,~o dti 


d 


dt. 


i+2k n +i — 1 


by the definition of 53* in equation ifT9l) , we have 
a +5+o 

= [55*,ad”; 1 P 0 ] 

oo oo . 

=(-i) n+1 E i +1 ” 2(Vl -E e n 5 


kj | U 


4/—2n 


5 


& 71+1 —1 


l=n ki,...,kn>l \j =1 

Ey=i*i=i 


0 t 


2 Z — | —2^7^,—)— ]_ —77/ 


' 72+1 


4-d" +i e e n 4 


4i—2(n+l) _ 


0 


i=n+l fci,—,fc n +i>l \j'=l 

E n- 

j =1 '“•> “ 

7772 1 


dt2i—n 


This proves the formula for ae?^ t Po 5 and we obtain 

OO 7 1 / 772 

p .=EE s £ (II s 


t » 


2 « 


0 


j=l m=l fci,...fc m >l \j=l 

E^Li 2 kj—i+m 

Then the lemma follows from equation dill) . 


dti 


i +1 


□ 


Note that, by equation f[9|) , we can obtain the operator P in Theorem [T] from the 
change of variable as: 

^ ~ ^ltfc=(2fc-l)H<?2fc+i ‘ 

Next, we define another linear map 02 from the space of power series of x and 
y to the space of second order differential operators with constant coefficients of 
variables ti by the following formula 


©2 (E = 

i,j > 1 


0 2 


i,j> 1 


a:?: 


dt dti 
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where ctij are constants. This map is an isomorphism when restricted to the 
subspaces 

d 2 

©2 : { I = «ii} —)• ( a './7)77]7 

M>1 i,j>l ' ' 

with the inverse map given by 

Using this isomorphism, we have the following lemma. 

Lemma 15 Define the series Q B (x,y ) to be 


Q B {x,y) = 


l-exp(®(i) + ®(±)) 


x + y 

I v^ c 


1 _ pxn J V°° B 2k ( r 2k-l , 2fc—1 
1 eX P \2-,fc=l 2fc(2fc—1) + +2/ 


x + y 


Then 

exp(«Bi + ^Qo') = exp(<Bi) exp(^Q^) 

where 

Qt = e 2 « B (*,i/))|, 

Proof: Using the Zassenhaus formula, we can see that 

°° / i \n—1 

<+ = E 

ra=l 

Similar to the induction we have made in the proof of previous lemma, we can 
also derive the formula for the nested commutator adfiffi 1 . First, for simplicity, 

we will use the following convention for summations and products: (• • •) = 0 

and nf (' • •) = 1 f° r l > k. Using the notation (EHl) . we have, 


= E II B ^n 

ki,...,k n >l \m =1 


2(2fe m -l) 


n —1 


E (-D‘ + "E 


i,j> 0 

i+j=2k\—2 


1=0 


n — 1 


d 2 


1 ' dt i+j: l +f 2 ^km-l) dt j+j:l=l+2^m-l) 
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We prove this formula by induction on n. In fact, when n = 1, 


ad% t Q™ = Q^ = Y1 B * u2(2fc_1) E 


i,j> 0 
i-\-j=2k—2 


dtjdtj 


Assume that, for n > 2, the above formula for ad^ t l Q q holds. Since 


“ 3 dtj +2 k n+1 -i ’ dt a dt b dt a dt b+2 k n+1 -i dt a+2 k n+1 -idt b ’ 


we have 


ad\Q™ = [K u ad n ^QV 


E II B k m ^ {2km ~ l) 

fci,...,fc n +l>l \m=l 


E m)‘ 


i,j> 0 

i+jr=2fci—2 


n — 1\ cr 

1 ' ^+E^ 2 (2fcm-l)^+E”t 1 i+2 (2fc m -l) 


n y 1 

E / 77- — 1 

i-i 


1 7 dt- . sr^l+l (^i i\9t • .v^n+l 


e n a 

fci,...,fc n _|_i>l \m=l 


E <-i) i+n+1 E 

ij>0 /=0 

i+ji=2/ci—2 


b+E£d 2 (2fcm -1) 17 "i+Er = A 2 (2fc m -1) 


,2(2fc m —1) 


1 ' ^+E^= 2 (2^-l)^+E™t 1 !+2 (2fc m -l) 


Note that in the above equation, we have used the Pascal’s rule 


n\ (n — 1 \ (n — 1 

l) = (l- l) + ( l 


for 1 < l < n — 1. This shows that the above formula for the nested commutator 

io Trolirl TTorieQ 


ad^ t Q o is valid. Hence 


oo . / n 

QV = E Ei E IT ^E (2fcm_1) 


n=l \^=1 


E m) ,+1 E 


ij>0 

i+j=2k\— 2 


n — 1 


1 ' ^+E^2(2fem-l)^+E”= (+2 (2fem-l) 
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On the other hand, by the definition of Q B (x,y) in the lemma, we have 

00 1 / 1 1 \ n 1 

nl \ x y 

n= 1 v y 

(®(j) + *(;)) 


x + y 


x + y 


E 

n= 1 


n\ 


«(-) + »(-) 
x y 


n —1 


Since x 2fe 1 + y 2fe 1 = (x + y) ^i=o 2 (~ B) l x l y 2k 2 *, we can expand Q B (x,y) as: 


/ 


Q B (x, y) = 


\ 


X> fcl (-!r + w 

ki=l i }< j>0 

i+j=2/ci—2 


00 1 n— 1 

E^E 

n=l «=0 


l J x y 


\ 
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\*+i h„i 


ki=l i,j>0 

i+j=2fci—2 

00 1 n—1 

E s E E 

n=l fc 2 ,...,fc n >l 1=0 


n — 1 
Z 


/+1 


n ^ 


„2fc m -l 


\rn=2 


II 5 fcm?/ 2fcm 1 

Vtti=/+2 


00 . / n 

e^ e n*. 


n\ . 

n=l ki,...,k n >l \m= 1 
n—1 

E M) <+1 E 


ij>0 

7+ < 7=2Ali—2 


/=0 


n — 1 
l 


j*+Em = 2(® ra l)yi+Em=I+2(^ m 1) _ 


For the degree of u in QY, observe that 

Z+l n n 

i + ^ " (2A: m — 1) + j + ^ " (2A: m — 1) = —1 + ^ ] (2 k m — 1). 


771—2 


771=/+2 


771=1 


Then can be obtained from © 2 (Q B (x, y)) by the variable change t k —>• rt 2k 1 t k . 
This completes the proof of the lemma and Proposition [H 

□ 


Furthermore, we can write Q B (x , y) as the following expression, using the Stirling’s 
approximation for the gamma function (|5|) : 


Q B {x,y) 


1 I i+l 

] xxyye x y 1 1 I 

—— —r(-r(-—. 

x + y ‘I'Xyfxy x y x + y 


(45) 
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We will use the above expression later in our proofs. From the definition of 
Q B (x,y) in the previous lemma, we can see that its expansion is of the form 

Q B (x,y) = Y QijX l y J , (46) 

i,j> 0 


for some constants Q B with Q B = Q B t . Now we check that the following action is 
well-defined: 

exp (^QY) exp (Pt) ■ exp (F K (t)). 

The operator exp (Pt) performs a shift —>• t/ ; — Ck~iU 2k ^ 2 for k > 2, which 
means that to and H remains unchanged, and other tk will be added by a term of 
u with positive degree 2k — 2. Hence exp (Pt) ■ exp (Fx(t)) is a well-defined formal 
series. Since the degree of u in any term of the operator Q] v is always positive, the 
coefficient of a fixed term in exp(4Q)' / ) exp(Pt) ■ exp (Fx(t)) can only be a finite 
summation. Hence the above action is well-defined. 

We go back to equation ©• By Lemma M, Lemma m and the discussion 
above, we have concluded the relation (1201) presented at the beginning of this 
section. 

Remark: We note that equation (1201) and the expression of Q B (x,y) are similar 
to results on page 11 of Givental’s paper j9j. Givental’s results also describe a 
relation between Kontsevich-Witten tau function and a generating function for 
Hodge integrals with a parametrization which is different from that of the Hodge 
tau function considered in this paper. The power series which determines the 
second order differential operator in his paper was deduced from the theory of 
Frobenius manifolds. In our case, the series Q B (x,y ) arises naturally from the 
Zassenhaus formula. It is interesting to observe that these two power series are 
essentially the same up to change of signs. 


□ 


To prove Proposition [TJ we can use Zassenhaus formula to obtain a decompo¬ 
sition of the following form 


ex P( a mU m Lm ) = exp (Y a m u m X m )exp(]-Q + ). 


(47) 


m> 0 


m> 0 


To compute the operator Q + , we set 

= Y a m u m X m and T + = a m u m Y m , 

m >0 m >0 

where X m and Y m are defined in equation (HTH) . Since X + € 0 i, Y + € 03 , we have 


1 

2 


71=1 


(- I )"" 1 

n! 


ad r Y Y+ - 


(48) 
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Such a decomposition is entirely similar to the case of exp {Y2 m <o a -mU m L m ). 
To see this, first, we define g,, g 3 to be 


01 — ^ Y^ oiijqi+j 


d_ 

dqi 


otij are constants 


) 


i Y ^b ( io/ib 

are constants, 7 ab = 7 fea 

a,6> 1 

J 


/ / 


By a similar discussion at the be 0 _ 0 ...._., ___ 01)51 

and g x ® g 3 are all Lie algebras. 

Let S be the Lie algebra isomorphism from g, ® g 3 to gi ® g 3 (in variables 
q t ) defined by 

d ^ _ i + j n d 

' Wi+j 0 j — Qi o > 

dqi i dq i+j 

and 

8 2 

S( qiqj) = -ij a „ • 
dqidqj 

The fact that E is a Lie algebra isomorphism is due to Lemma Q33 in Appendix 
A.2. 

For m < 0, we also write L m = X m + Y m where X m is defined as in equa¬ 
tion m and 

1 


Ym — 2 'y y QaQb- 


Set 


a+b=—m 

~ m Xm and Y~ = Y a -m,u~ m Y m . 


X~ = Y a -mU~""X m 

m <0 m<0 

We can see that 

S(-X~) = X + and E(-y) = Y+. 

Also, by the Zassenhaus formula, we have 

exp(— Y a -mU~ m L m ) = exp(-A _ )exp(-^Q') 

m< 0 

where 




°° ( i \n —1 

= E L ^, —uT-xU-Y-). 

n =1 

Since E preserves the commutator relation, we can easily deduce that 


1 {-Q~)=Q H 


(49) 
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Inspired by Lemma 4.5 in [TT], we introduce a method to calculate the quadratic 
form Q~ explicitly, and then obtain the form of Q + using the isomorphism E. 

Suppose {V irL } is a set of differential operators indexed by negative integers 
satisfying 

[bm.) ^n] = {pi ^)bm+n- 

The first thing we want to do is to differentiate the operator exp(A(it)) by u, 
where 

A{u) = ^2 a mU~ m V m , 

m< 0 

and a m are constants. Then we have 


Lemma 16 For a fixed sequence of numbers {a m }, there exists a unique sequence 
{dm}, such that 


®Mu) 

8u 


{^2d m u- m - l V m )e A M. 

m< 0 


Here we are dealing with differentiation of an operator-valued function. For more 
details about this and a proof of the above lemma, we refer the readers to Appendix 
A.3. Observe that, for negative integers m, the sets of operators {L m }, {X m } and 
{z 1 ~ m d/dz} all satisfy the same commutator relation as {V m } defined before. 
Hence they all satisfy Lemma [T6l 
Let 


E{u) = exp(^ a m u m L m ) 

m< 0 

= exp(^ a m u~ m X m )exp{^Q L {u)), 

m< 0 


where Q L (u ) is a quadratic form. Then, 

^ E(u ) = | ^ ex P(^ a m u~ m X m ) j exp (^Q L (u)) 

\ m< 0 / 

+ exp(^ a m u~ m X m ) f-^-exp{^Q L {u)) \ . (50) 

m< 0 ' ' 

Consider the first term on the right hand side of equation (l50l) . Since {X m } satisfy 
Lemma flHl we have 



exp(-Q L (ii)) 



E{u). 
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For the second term on the right hand side of equation (1501) . its action on an 
arbitrary function G in variables qi is 


exp(-Q L (ri)) ) )■ • G 


| exp(^ a m u m X m ) 

\ m< 0 

= exp(^ a m u~ m X m ) • | f^-exp {^Q L {u)) \ G 


= exp(J] a m u m X m ) ■ | \Q L( y u )) (exp(^Q L (u))G^ | 

= | exp(^ a m u~ m X m ) • -^\Q L ( U ) | | ex p(^ a m u~ m X m ) ■ fexp(^Q L ( u))g\ 

l m<0 J l m<0 ^ ' 

= I exp(^ a m u~ m X m ) • | £(«) • G. 

I m<0 U J 

Recall that we use “ • ” to represent the action of differential operators on func¬ 
tions. The third equality of the above equation follows from the fact that X m is 
a derivation. Since G is arbitrary, we have the following equation for differential 
operators 

exp(^ a m u~ m X m ) (^exp CjQ L {u)) 


m< 0 


l exp(y~] a m u m X m ) 

\ m< 0 


d_l 

du2 


Q l {u)\E(u). 


Hence, 

d__ 

du 


E(u) = ( ^2 d m u m 1 X m + exp(^ a m u m X. 

\m <0 m <0 


m> ' du2 QL<<a) ) E ^' 


On the other hand, applying Lemma fl6l on the case V m = L m , we have 


^-E(u) = ^ d m^ m ^Lm)E(u). 


m< 0 


Therefore, since E(u) is invertible, 


= exp(— ^2 a mU m X m ) • ^2 d 


-m-ly 

, m u, 1 n 


(51) 


m< 0 


m< 0 


This gives us a way to compute Q L (u ) when the coefficients {a m }, which also 
determines {d m } by Lemma fl6l are given. Note that by the definition of Y m , YLi 
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is defined to be 0. And we can see from the above equation that the power of u 
for the term qiqj in Q L (u ) is always i + j. Usually, an easy way to obtain { d m } 
is to apply Lemma [T6l on operators {z 1 ~ m d/dz}, when we know the explicit form 
of the function x(u, z ) defined as 


x(u,z) = exp(^ 


Again, by Lemma fl6l 
d 


B~a eXp{ ^ 


m< 0 


—m _,1 —m 


a m u- m z x ~ m ^)-z. 


a m u z 


m< 0 


d_, 
dz' 


n.-m -1 1 -m 

HI a * 


d_ 

dz 


'exp( J] 


—m _1 —m 


a m u z 


d_ 

dz 


(52) 


m< 0 m<0 

We use the operators on both sides of the above equation acting on z. Since 


d_ 

du 


exp 


m< 0 


a m u- m z 1 - m ^-)-z = 

dz du 


<P (2 


—m 1 —m 


exp( > a m u z 

m< 0 




z 


we have 


d 


— x(u,z) = (J2 d r 


—m— 1 —m 


9 


m< 0 


* )&*(“’*)■ 


(53) 


So the coefficients {d m } can be easily obtained by comparing the partial derivatives 
of x(u, z) with respect to u and z. Now, for the sequence of numbers {a_ m } 
determined in Lemma fTTil m < 0, let a m = —a_ m , then 


Q L (u) = -Q . 


Now, to determine Q , it is convenient to consider the power series Q(x, y) defined 
by the following function 


Q(x,y) = log ( 


1 


1 )(. xy 


h(x ) h{y) y — x 


(54) 


where h is the function defined in Corollary 1121 The series can be obtained in the 
following way. Using the series 1/h in Corollary 1121 we have 


1 


1 00 7—2 7—2 

=xy^-iy-'ibi y 


h(x) h{y) y-x 


i= 1 
oo 


y-x 


i—2 




i= 3 


m =1 
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Since the series l) l ibi Y^ m =i x m y l 1 m has no constant term, the loga¬ 

rithm of the above series is well-defined and we have 


Q(x,y ) = 


(- 1 ) 


n —1 


i -2 


n=l 


n 


Et- 1 )^ E 


a=3 


m=l 


Let the series Q(x, y) be of the form 


Q(x,y) = E ® 
i,j> 1 




(55) 


for some constants Qij. From the definition (IMP we can see that Q(x, y) = Q(y, x). 
Hence Qij = Qji. Next, similar to the map 02 , we define the linear map ©3 from 
the space of power series of x and y to the space of quadratic functions with 
constant coefficients in variables qi by the following formula 


©3(E = E "b'M.b 

i,j> 1 i,j> 1 

where are constants. This map is an isomorphism when restricted to the 
subspaces 


©3 


: { Y a ij xl y J 

i,j> 1 


a. 


V 


— CVji} 


(E a 

i,j> 1 


V 


QiQj } 


with the inverse map given by 

© 3 1 ('Mj) = ^(®V + ®V). 


Now we prove that 


Lemma 17 

Q~ = ®3(Q(x,y))\ qk ^ u k qk . 

Proof. After setting a m = — a_ m , the function x(u,z) defined in equation (1521) 
becomes r] = 77 ( 11 , 3 ) defined in Corollary [T3j Since 

f) n 7 2 1 

—n(u z) = -- -i_-_- 

du ^ u u(l + uz) 2 r] 


and 


b (u ' z) 


z 1 

(1 + -uz ) 2 r]' 
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we have 


lr ,( “, 2 > = (-- 


2 (1 + uz) 2 2 \ d 


u 

OO 


uz ■" 


E + l )”- 1 


5 


\n=2 


(n + l)n(n — 1) 


U Zn )Wz r,M 


From equation 


we can see that, in this case, 


d— n+l — ( 1) 


71—1 


for n > 2. Let 


(n + l)n(n — 1) 


Q = Qij ul+3 9iQj 

i,j> 1 


for some constants Q- with (+• = Note that as explained in the remark after 
equation (15T1) . the power of u before q^j should be i + j. Hence 


ld_ 

2 du 




i,j> 1 


On the other hand, by equation (15111 . 


l~L Q = e M^2 a ~ mU ~ 

m< 0 


1 X m )• Y, d -n+lU n - 2 Y. 


—n+l 


\n =2 


We want to find a series T(x,y), such that 


n—2 




n=3 


i=1 


which gives us 


(56) 


n—2 


2/)) — ^ ^ n+l ^ ' QiQn—l—i — /* ' n+1^- 

71=2 


—n+l- 


n=3 i=l 

This series is nothing but the following one 


x ,1+?L 


2/ ,1+s, 


1 1 3 


T(x,y) = -(-)"log(l + 2 /)-(-) log(l + s) 

x — y y x — y x x y 2 
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One way to see this is that 


OO 

x — y s—' 1 


n —2 


xy *■—' 2 

y n =3 


Ei - 1 )" -1 


EE-+iE*V* -1-i 


2=1 

2 


77=3 

OO 


E )- 1 )" -1 


(n + l)n(n — 1) 

1 2 1 


77=3 


2 

-+ 


n + 1 n n — 1 


(x n - 2 -y n - 2) 

( x -2 _ y n - 2 ) 


=-o (log(l + x) - x + — - —)- 2 ( 1 °s ( 1 + x ) - x + Tr)-Og(l + - x) 


x 


+ -4 ( lo §( 1 + y)-2/ + 4 _ 4) + 4( 1 °s( 1 + y) ~ v + 4) + — 0og( 1 + v)~ y) 

yo 2 6 y z 


1,1 + ^2 


=-(-) log(l + y) - -( 


1 .1 + X, 


1 


2 ' y 
1 3 3 


y y 

x-y 


x x 


log( 1 +x)-^ + ^-- + - 


xy 


T(x, y). 


Note that T(x,y) = T(y,x). Similar to the isomorphism ©i defined in Section 
3.1, we define 

OO OO OO 

0 i otiX 1 ) = ©i(^ a t y l ) = 22 a iQi 
2=1 2=1 2=1 

for constants a%. In fact, it is easy to see that, for two series fi(x), f 2 (y) 


®3{fi(x)f2{y)) = &i(fi{x))@i{f 2 {y))- 

On the other hand, since h l (z) = e®* ■ z l by the definition of h in Corollary 
fl2l by equation we have 

@i(h l (x))@i(h j (y)) = exp(^ a- m X m ) ■ (q^qj) 

m< 0 

for all i,j > 1. This implies that 

( OO 77—2 \ 

E 5 < 1 -n+lE ,l '< ;C )' 1 ” -1- ‘(!') 

27=1 Z 2=1 / 

( OO ^ 77—2 \ 

^ ^ 2 ^—77+1 ^ ^ QiQn—l—i | • 
27=3 2=1 / 

Then, by equation (|56l) . and setting u = 1, we have 

0 3 (T(h(x),%))) =~ 22 (* + = “© 3 ( 22 + j)% xl y 3 )- 

i,j > 1 i,j> 1 
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The series h has no constant term. Hence the composition T(h(x),h(y )) is well- 
defined. Since both sides of the above equation are symmetric with respect to x 
and y, we obtain 

T(h(x), h{y)) = + j)Q~ jX l y 2 . (57) 

i,j> i 


To simplify the expression for T(h(x), h(y)), we take derivative of both sides 
of equation ([551) with respect to 2 and replacing e~^ z by the left hand of the 
equation (1351) . We then obtain the following formula for the derivative of h(z ): 


d h _ z(l + h) 2 
dz h 

Also, after taking the logarithm on equation (1351) . we can obtain an expression for 
log(l + h), 

106(1 +h) = r 2 + TTh- 

Therefore, 


(1 + h ) 2 
h 2 


log(l + h) 


Using the above equation, we obtain 


lid h 1 _ 

2 ft 2 dJ + ft + L 


T(h(x),h(y)) 


h{x ) 

h(x ) - h(y) 


y dh (y ) , _J_ 

2 / 1 ( 2 /) d y h(y) 


_ h(y) f x d h{x) 1 \_1_1_3 

h{x) — h(y) \2h(x) dx’ h{x) ) h(x) h(y) 2’ 

h[x) y d h(y) h(y) x d h(x) 1 

h{x) — h(y) 2/i(?/) d y h(x) — h(y) 2 h(x) dx 2 

h{x)h{y) f y d h(y) _1_ _ x dh(x) 1 \ 

2(h(x) — h(y )) \h 2 (y) dy h(y) h 2 (x) dx h(x)J 


By equation (1571) . we have 


Qij^ + j) xl y ] = 

i,j> 1 



1 N -1 ( d(x//t(x)) 

Hy) V dx 


d {y/h(y)) \ 
d y ) 


We replace x and y by ux and uy respectively in the above equation. By the chain 
rule, we have 

+ jV+^x V = ( 

i,j> 1 



U 


U 


X -1 5 


h(ux ) h(uy) du \h(ux) h(uy) 


u 


u 


(58) 
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By the formula of 1/h in Corollary 1121 we have 


u 


1 1 


h(ux) h{uy ) x y 


+ £(-l ) n - l nb n u n -\x n ~ 2 - y n ~ 2 ) 


n =3 


n —2 


i 1 - - b 1 +i . 

x y \ ri=3 t=i / 


Let 


K(x, y, u) = (---) 


-l 


u 


xy / tt 


x y \h(ux) h(uy) J y — x \h(ux) h(uy) J 


u 


Then 

oo n— 2 

il (x, y,«) = 1 + ^ xy- 1 -’ 

n=3 i=l 

and logiL(x,r/,n) is well-defined. Equation (15S1) can be written as 

Qij^ + j)u l+3 ~ lxl y 3 = K(x,y,u)~ 1 -^I\{x,y,u) = J^log K(x,y,u). 

i,j> 1 

This shows that both series Qi 3 u l+3 x l y ] and log K(x,y,u) have the same 

derivatives with respect to u and both of them are 0 at u = 0. Therefore we have 

Y Q7j ul+3xl y J = log K{x,y,u) = Q(ux, uy). 

i,j> 1 


Applying ©3 on both sides of this equation, we obtain the lemma. 

□ 


Remark: Until now, we have introduced two series Q B (x,y ) and Q(x,y). In this 
paper, we use the formal power series to describe the coefficients of our differential 
operators, and we only need to capture the values of the coefficients in the series. 
The values of x and y are not important to us. 

□ 


Using the relation P9l) . we can deduce that 
Corollary 18 The operator Q + in equation (1471) is given by 


Q + 


Y Qij ul+Ji J 

i,j> 1 


d 2 

dqidqj ' 


This completes the proof of Proposition [7) 
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3.3 Proof of Proposition |8] 

To prove this proposition, it is equivalent to show that 

Qfj = (2* + l)!!(2j + l)!!Q 2 i+ 1 , 2 j+ 1 - 

where Q B are the coefficients of the series Q B (x,y) defined in equation (TTOD and 
Qij are the coefficients of the series Q(x,y) defined in equations (|5iD and ([55]). 
i.j > 0. Let 


d 

A(x,y) = x—Q(x, y) 
oy 

d 


=V og|( 


1 .)(. xy 


h(x) h(y) K y — x ) 


Since h = w 1 — 1 by equation ([Mil . we have 

w' (y) 


A (a;, y) = x- 


x xw (y) x 

+ r- 44 + -• 


w(y) — w(x) x — y 1 — w(y) y 
Using the Gaussian integral 

y/2n J- 

We consider the following integral transform 


(59) 


r oo 2 

z 2k+2 e -^ dz = (2k + 1)!!, k > 0, 


I(x,y) = 


2vr 


_ z -ps 

A (xt,ys)e 2 dfds. 


On one hand, from equation (1551) . we can see that the expansion of A (x,y) is of 
the form 

A(x,y) = jQijX^y 3 ^ 1 . 


i,j> 1 


Hence the expansion of I(x,y ) is 


I (x,y) = - V Q2a+i,2b+i(2a + l)!!( 2 b + l)\\x 2a+1 y 2b+1 . 
11 

u a,b >0 


(60) 


On the other hand, for the last two terms on the right hand side of equation (15911 . 

/(f 


xtw (ys) xt. t 2 

-44 + - e"Td t = 0 

- w{ys) ys 


since the integrant is skew symmetric with respect to t. Therefore we only need 
to consider the integral transform of the first two terms of equation (15911 . namely 


I (x,y) = ^(h + I 2 ), 
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where 


h = / / xtw ( ys) 


1 


w(ys ) — w{xt) 


_ tf+£_ 

e 2 dtds, 


/'2 = 


xt _£+s£ 

-e 2 dids. 


xt — ys 


We first calculate the double integral I\. Let w\ = w(xt),w 2 = w(ys). From now 
on we consider w\ and W 2 as functions in variable t and s respectively, and treat 
x and y as non-zero constants. By the relation 

x = w (x)(l — w~ 1 (x)) 

in equation (1141) . we can transform I\ into the following form 


h = 


1 dwidw 2( 1 . 1 _ * 2+ « 2 ,,, 

-(1-)-e 2 dtds. 


xy dt ds w± w 2 — w\ 


Observe that the two series v and w introduced in Section 2.3 have the relation 
w(—z) = v(z). We extend the domain of w to R by defining w(z) = v(—z) for 
z < 0. Then the range of w is (0,oo). Also w(xt ) satisfies 


w(xt)e 


1 -w(xt) _ / g-V 


by equation (USD- Hence 


_9 _9 _9 

e~T = wf e x ~ W1X . 


The integral I\ can be written as 


x 2 +y 2 roo /*oo 1 1 

h = - - / / (1-)- wf V *e~ wix *- W2V ~dw\dw 2 - 

Jo Jo 


„-2 „-2 2 .„„„,-2 


xy Jo Jo 


W 2 — w 1 


Consider the change of variables (wi,w 2 ) —>• (z,u),w\ = x 2 zu,w 2 = y 2 z(l — u). 
Since w\,w 2 > 0, and 


-2 -2 w i x 

z = w\x + w 2 y , u = 


-2 


w\x 2 + w 2 y 2 ' 


we have 2 E (0, 00 ) and u E (0,1). Then 


, 1 = LA!L (l2r+vr+I 

xy 


. (l-^—)^—^ Y ^-z x ~ 2+y ~ 2 e- z dzdu. 

/0 Jo + y*)u 
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We solve the double integral in the above equation by two steps. First, we integrate 
by z. Using formula ©, we have 


u x ~ (1 - uf 

y-2 _ ( x 2 J o 

u x ~ 2 ~ l {i-u)y 


roo 

/ (1 

J o 


X 2 ZU 


)z x 2+y ~ 2 e~ z dz 


u 


y 2 - (: x 2 + y 2 )u \ J 0 
u x ~ 2 -\i - u y~ 2 


z x 2+y 2 e~ z dz- 


—2 


r 

Jo 


x- 2 +y ~ 2 


y 2 — ( x 2 + y 2 )u 
{x~ 2 + y~ 2 )u — x~ 2 


(■ uT(x 2 + y 2 + l)-x 2 T(x 2 + y 2 )) 


r (x- 2 +y~ 2 )u x “ A (l -u) 


y 2 — ( x 2 + y 2 )u 
= - x - 2 y _ 2 r(x -2 + y~ 2 )u x ~ 2 -\l - u) y ~ 2 . 



Secondly, we integrate by u. Using formula (fTElh we have 


/ 


u 


L (1 — u) y ~du = B(x 2 ,y 2 + 1). 


Hence, 


h = - eX +y (x 2 r 2 (y 2 r 2 nx- 2 +y- 2 )B(x- 2 , y - 2 + 1 ) 

xy 

= 2 )„-» r( ^ 2)r( »- 2) , 
y x z + y z 

where the second step above is from equation (fTTl) . For the second integral I 2 , we 
consider the change of variables u = t — ^ s,z = s. Then /o can be transformed 
into 


h 


u+^z 

-—e 2 e 2 dudz 

u 

(u+%z) 2 y f f _i 

e 2 e 2 dudz H— zu e 2 e 2 dudz. 

x J J R 2 


Next, we will use the following two formulas of the Gaussian integral to calculate 
the above integral I 2 : For a > 0, 



[tt &£. 
\ —e 4a 

V a 


and 



Xe -ax 2 -bx dx 


\pKb _3 &£. 

- a 2 g 4a ' 

2 
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First, we compute 


(«+ g *) 2 _ z £ 
e 2 e 2 dudz = 


/ OO 2 ^2,2 /*00 

e-¥(^) r 

-OO J — 


-W-lzu dudz 


e 2 X 


. 


00 z 2 

e 2 " dz 

— OO 


=2vr, 


and 


_Cu+| f F _£ ^ 

zii e 2 e 2 ditdz = 


Z " 00 —i —— f 

/ u e 2 
J —oo J — 


2 /-OO *2,2 2 


£l±3CZl_3£„ 


ze a: 2 2 x a dzdu 


= - ^( x2 + ;,2 rl 


e :r2 +!/ 2 2 dtt 


= -2vr 


x x x 

xy 

x 2 + y 2 


Then we have 


Finally, we obtain 


y xy x 

I 2 = 2n — 2n — — -- = 27T- 


x x 2 + y 2 x 2 T y 2 


—I( x i y) + I2) 

x 2ir x 


xy 


1 (x 2 Y 2 (v 2 ) y 2 e x 2+y 2 1 1 

- j V, 2 - r (—)r(—) 

x z + y 2 x z y 2 


x 2 + y 2 27t 
=xyQ B (x 2 ,y 2 ), 


where the last step above is by equation (H5l) . Using equations (1601) and (Rbl) . we 
can deduce that 

(2a + 1)!! (26 + l)!!<32a+l,2b+l = Qab- 

This completes the proof of Proposition 0 and thus also completes the proof of 
Theorem [I] 


□ 


4 Proof of Corollary 2 

We prove Corollary [2] in this section. Our following discussion is quite similar 
to the one in [2], except that it is based on Theorem |T] we have proved, and we 
also use the Zassenhaus formula, which is more clear in terms of the computation 
here. Later we can see that, since equation (1631) holds for any set of coefficients 
{l m }, the expression similar to equation d2|) that connects the two tau-functions 
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using exponential of Virasoro operators and first order differential operators is not 
unique. 

The Virasoro constraints for the Kontsevich-Witten tau-function in variables 
t are (cf. [13]): 


(2fc + l)!! _d_ 

m L, (2k-2m-l)\\ tk - m dt k + 2 

k>m k-\-L=m—l 


\ E (2* + l)!!(2* + l)!! 


d 2 


dt k dti 


- (2m + 3)!! 


d 


tn 


dt 


m-\-1 


+ ~^5 m -1 + (' m — 0:1; 2 ,...) (61) 


such that 

L m • exp(F K (t)) = 0, 

for m > — 1. If we remove the terms in L 2 m involving the even variables q 2 k , 
then it coincides with the sum of the first two terms of the operator L m after the 
variable change = (2k — l)!!<72fc+i- Now, we choose the Virasoro constraints L m 
with m > 1. Since exp (Fx(q)) has no even variables q 2 k, we can deduce that, 

V 2m ■ exp (F K (q)) = 0, 


where 

~ d 

^2 m = L 2 m (2m + 3)tt . 

Oq2m+3 

Hence, for any sequence of numbers {l m }, we have 

oo 

exp( E lm V 2 m) ■ exp (F K (q)) = exp (F K (q)). 

m= 1 

Using the Zassenhaus formula, we have 


exp( E lm V 2 m) = exp(E l mL 2 m) exp(- E h k 


d 


m= 1 


m= 1 


k =1 


dq 2 k +3 


), 


for a set of numbers {bk} with 
bi = 5Zi; 


n> 2 


n—1 


bk — (2k + 3 )lk + (2k + 3) E s E Imi (3 + 2m \)... l 77lri (3 + 2 E m j). 

mj >1 j =1 


m j= k 


(62) 


The above formula can be obtained using the methods introduced in the proof of 
Lemma 1141 The induction is entirely similar. So here we will skip the details, and 
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only give the formula for the nested commutator 


U/U/r-\00 7 T 

l^m= 1 L m^2r 


t ( — lm(2m + 3) 


5 


ra=l 


dq2m+3 


) 


yi ^mi(3 + 2mi)... lm n (,3 + 2y>,) 
fc=l 77ty>l j = l 

E”=i"*i=fe 


5 


dq2k+3 


We can also see from equation (1621) that {6fc} uniquely determine the coefficients 
j7 m } and mce versa, using the same deduction in the proof of Lemma [121 And 
replacing l m with u 2m l m , we have 


exp(^ l m u 2m V 2 m) 

m= 1 


oo oo 

exp( ^2 lmU 2m L 2m ) exp(- ^2b k u 2k 
m=l fc=l 


d 


dq 2 k+3 


Hence, we can conclude that 


f) 

ex P( - X] ^ ku2k dq 2k+3 ^' exp '( Fk ) 


= exp(— ^ l m u 2m L 2m ) ■ exp (.F^ (<?)). (63) 

m= 1 


In our case, if we set b k = b 2k+ i, then the operator exp(P) in Theorem |T| can be 
replaced by the exponential of Virasoro operators at the right hand side of the 
above equation with the corresponding coefficients {lm}- Another way to represent 
the relationship between the two sets of coefficients {b 2k +i} and {l m }, like we have 
done before, is to use a power series. Let 


9(z) = exp(— ]T l m z l ~ 2m 

m= 1 



1—2 m 


d 


—3 


Then, by equation 

OO 

0~ 3 (z) =exp(- ^2 q z j ~ 

m= 1 


n —1 


= 2 : 


n\ 

n> 1 


• lrn n (3 + 2 ^ ^ uuj 
i=i 


.-W 3 2 E? = i m 3 


And by equation (1U2D . we have 


0 - 3 (z) = 3^ 

k =0 


fofc+l „-2fc-3 

2fc + 3 
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which gives us 


0{z) 



^2k+l 2fc—3 

2fc + 3" 


Note that if we transform the above equation into 


1 

3 


1 _ = hk+l —2k—3 

3 6 3 {z) 3z 3 ^2k + 3 

k=1 


this expression is very similar to equation (2.146) in [2|. The first few terms of 
9(z) are 


fln l-i 67 

9(z) = z - z - 

v ’ 180 453600 


z 3 + ... 


Furthermore, the following two operators can be combined into one using the 
Baker-Campbell-Hausdorff formula: 


oo oo oo 

exp(^ a m u m L m ) exp(— ^ l m u 2m L 2m ) = exp(^ e m u m L m ). 

m= 1 m= 1 m= 1 


Since {u m L m } and {z 1 m d/dz } satisfy the same commutator relation as V m , that 
is, \y m , 14,] = (m — n)V m+n , it is easy to see that 

OO OO n OO 

exp(^ a m z x ~ m — )exp(— ^ = expe m z 1 - m —). 

ra=l ra=l m=l 


Using #(z), we can obtain the coefficients {e m } from the following series 


d 


d 


exp(^ emZ 1 m —) ■ z = 6 ( exp(^ cimZ 1 m —) • z ) = 0(/(z)), (64) 


m=l 


ra=l 


<9z" 


where /(z) is defined in equation (12811 . This proves Corollary [2j The first few 
terms of 9(f(z)) are 


«(/«)=.+§-s*- 1 +ik a+ -- 

The computation here is very straightforward, since both series 9 and / are given 
explicitly. We have computed the first three coefficients e m by hand, they coincide 
with the data listed in Alexandrov’s paper [T] (with the opposite sign, since our 
operators act on the Kontsevich-Witten tau-function, not Hodge tau-function): 


m 

1 

2 

3 

Cm 

2 

_3_ 

_PT 

2 

PT 
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5 Further remarks 


In the paper [2], Alexandrov established the following relation between exp q )) 
and exp(F^-(q)): 

exp (F K (q)) = c(u)G + ■ exp (F H (u, q)), (65) 

where c{u) is an unknown Taylor series which is conjectured to be one, 

OO 

G+ = /H L °exp (J2a k r k L k )^ L \ 
k=1 

where the coefficients {a k } are determined by a series /+, and u = /3 1 / 3 . More 
precisely, {a^} are determined by the relation 


d 


f+=exp(^2a k z 1 k — , 


k=L 


and /+ is given as a solution of the following equation 


f+( z ) 

1 + f+( z ) 


exp( 


f+( z ) 

1 + f+( z ) 


E exp (-E), 


where 


E = 1 + 


1 + f+( z ) 


) 2 



The operator j3z L ° is treated as the exponential of the degree operator | log(/3)Lo- 
The action is 

exp(^ log {/3)L 0 ) ■ q k = /3^ k q k . 

The use of parameter (3 in the conjectural equation (1651) is due to the relation 
between the Hodge and Hurwitz tau-function HU- 

At this moment it is not clear what is the relation between / + and 


OO 

ex p(- ^2 

m= 1 




1—m 


d_ 

dz 


)- Z . 


Based on our known results, it is possible to derive an explicit equation that has 
6(f(z)) as its solution. But we believe that the formula ([2]) may contain more 
information due to the appearance of numbers C k = (2 k + l)!!&2fc+i- If we switch 
to variables ti in P, the generating function after the shift 

exp(P) • F K (t) = F(t 0 ,t 1 ,t 2 - Ciu 2 ,t 3 - C2U 4 ,...), 


which also satisfies the KdV hierarchy, looks particularly interesting. We think 
that there may exist some special connection between this shift and the one derived 
in [2l]. 
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Appendix 

A.l Special cases of Zassenhaus formula 

The Zassenhaus formula is not as well known as the BCH formula and it needs 
more caution when applying Zassenhaus formula. In this appendix, we will verify 
formulas (1401) and (1411) . which was deduced from Zassenhaus formula before, using 
the following form of the BCH formula (cf.|19|): 


^e«) = X + Q -(f X ± 

V 0 n ( n + 1) J 


( 66 ) 


Note that for the summation inside the integral above, adx and adQ are treated 
as two non-commutative variables. Hence the integral will give us a formal power 
series in adx and adQ. Suppose X € fli and Q E 02 ( 01 ' 03 ). Then 

e adx Q E 02 (or 0 3 ), 


and, for any Q' E 0 2 (or 0 3 ), 


g wq = Q 


This implies that, for n > 0, 

( e ad A - e t a d Q ) n Q = ( g ad xjUQ 

Hence, in this case, the formula (1661) can be simplified into 

(1 - e ad *)A 


log(e x e«) =X + Q-[Y j 


—( n(n + 1 ) 


Q 


\n 
( 00 


^ + q-lE ( 1 -r ) ’ 1 -E ( 1 ~ e - ) " i q 


\n=l 


= X + Q-( — log(e adx ) + 


n= 1 

1 


n + 1 


log(e adx ) + 1 ) Q 


1 - e adx 


= X + Q- 
= X - 

= x + 


1 - e ad * 
adx 


adx + 1 ) Q 


1 - e ad * 

adx 


1 - e -adx 


adxj Q 
Q. 


Now we consider 


00 ( _ 1 \n —1 

Q = Y^ -- - - ad n v ~ l Y 


n= 1 

1 - e~ adx 


adx 


n! x 
lx\ 

Y, 
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for Y € 02 (or 03 ). Then 

,, X 0\ ( & dx 1 — e -adx 

log(e x e° = A + ( ----- 

\ 1 — e adx adx 

= X + Y. 

This shows that the two formulas (14011 and (1411) are valid. 


Y 


A.2 The Lie algebra isomorphism S 

Lemma 19 Let S be a linear map from 0 ^ ©03 to 0 i ©03 (in variables qt) de¬ 
fined by 

d i + j d 

1 — \Qi-\-j o / ■ Qi o 5 

oqi 1 dq i+j 


and 


S( qiqj) = ~ij 


d 2 


dqidqj 


Then E is a Lie algebra isomorphism. 


Proof. It is easy to see that the map E is a bijection by the definition. To check 
that E preserves the commutator relation, for example, we can check the following: 
Since 

dd d 

[Qi+j ~n ) Qi J = Qi+j ' 


I Q III 9 J 

oqi oqi-k 

d 


dqi-K 


[Qi+j-Q^iMk] = Qi+jQk, for i fc; 
r d Qi 

[Qi+j ^7 > Qi J = ^qi+jqi, 


by the definition of E, 




i + j d 

-Qi 


d 


HI n ) , Qi — k r\ 

i oqi+j I ~ k oqi 

i + j d 
■ i Qi — k o 
i - k oqi+j 

d 


oqi-k 
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and 



—“( Qi+jQk)'i 



— ^{^Qi+jQi) ■ 

The above calculations shows that H is an isomorphism. 

□ 

A.3 Differentiation of an operator-valued function 

The differentiation of an operator-valued function falls into the field of non- 
commutative calculus of operator derivatives. Derivatives of such functions have 
been discussed in the paper [23], where operators considered are those acting on 
Banach spaces. Since we are not working with Banach spaces, we will not directly 
use formulas in |24j . In our case, we treat A = A(u ) as a formal power series in 
u. We denote the formal derivative of A(u) to be A , 



m< 0 


First we show that 



Similar results can be found in [23]. Since our setting is different from that in [23] , 
here we provide a proof of the above formula for completeness . 

A straightforward computation shows that 



n =1 



n= 1 k =1 
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Since adjpA' = 2A{ad A A^)—ad 2 A X, and [A, ad A ] = 0, we can obtain, by induction, 

ad Ak A’ = ( A k - (A - ad A ) k ^j A'. 

Using the following formula for the rising sum of binomial coefficients 



for any i > 0, n > i, we have 


r\ oo -t n 

^ = E s E- 4 ‘- 1 a-«<M”-'A 


, n. 

n =1 k =1 

oo n n—k 




n\ 

n= 1 /c=l 2=0 

oo . n—1 


E^E^'E 


n\ 

n= 1 2=0 
oo n—1 


k=1 


n — k 


n — k 


(- 1+4 A 


(-1+4 A 


E s E ' 4 ”" 1 

ti =1 i 
oo n—1 

EE 


77 


1 




= e 


71=1 1=0 

OO OO 

E E 

2=0 71 = 2+1 
OO 

A 


(n — 1 — *)! (z + 1)! 

1 (- 1 ) 


A n 


(n —l — i)\ 

<E++A 


(* + !)! 


a+A 


Furthermore, using the formula 


e adA M = e A Me~ A , 


and the equality 


for any m> 1, we have 


E 

i,j> o 

2+J =771 


dd 


(l - i) m 
ml 


= 0 


_ 9 _ 

du 



(-+ 
d(j + 1)! 



e 


A 



e 


A 
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Proof of Lemma 1161 Since {V m } satisfy the relation 

[V m , V n ] = (m- n)V m+n , 


□ 


we can deduce that 



for some constants {d m }. Note that adP^ l X consists of terms u~ m ^ 1 V m with m < 
—n. Hence X^ m <o dmU~ m ~ l V m is well-defined and {d m } are uniquely determined 
by {a m } and the bracket relation of {V m }. 


□ 


References 

[1] A. Alexandrov, From Hurwitz Numbers to Kontsevich- Witten tau- 
function: A Connection by Virasoro Operators , Letters in Mathemati¬ 
cal Physics, January 2014, Volume 104, Issue 1, pp 75-87, Stable URL: 
http://dx.doi.org/10.1007/sll005-013-0655-0 

[2] A. Alexandrov, Enumerative geometry, tau-functions and Heisenberg- 

Virasoro algebra, Commun. Math. Phys. 338 (2015) 195-249, 

|arXiv:1404.3402vl [hep-th]] 

[3] F. Casasa, A. Muruab, M. Nadinic, Efficient computation of the Zassenhaus 
formula, Computer Physics Communications, Volume 183, Issue 11, Novem¬ 
ber 2012, Pages 2386-C2391 [arXiv:1204.0389v2 [math-ph]] 

[4] R. M. Corless, D. J. Jeffrey, and D. E. Knuth, A sequence of series for the 
Lambert W function . In Wolfgang W. Kiichlin, editor, Proceedings of ISSAC 
’97, Maui, pages 197-204, 1997 

[5] E. Dynkin, Calculation of the coefficients in the Campbell-Hausdorff formula, 
Doklady Akademii Nauk SSSR (in Russian) (1947) 57: 323-326 

[6] L. Comtet, Advanced Combinatorics: The Art of Finite and Infinite Expan¬ 
sions, rev. enl. ed. Dordrecht, Netherlands: Reidel, p. 267, 1974. 

[7] E. Date, M. Jirnbo, M.Kashiwara and T.Miwa, Transformation groups for 
soliton equations-Euclidean Lie algebras and reduction of the KP hierarchy, 
(pp. 635-653). RIMS preprint 362 (1981). 


51 



[8] C. Faber and R. Pandharipande, Hodge integrals and Gromov-Witten theory, 
Invent. Math. 139 (2000), no. 1, 173-199, |arXiv:math/9810173[math.AG]] 

[9] A. Givental, Semisimple Frobenius structures at higher genus, Int. Math. Res. 
Notices 23 (2001) 1265-1286 |arXiv:math/0008067 [math.AG]]. 

[10] A. Givental, Gromov-Witten invariants and quantization of quadratic Hamil¬ 
tonians, Mosc. Math. J. 1 (2001), 551-568 arXiv:math/0108100 [math.AG]]. 

[11] I. P. Goulden and D. M. Jackson, Transitive factorizations into transpositions 
and holomorphic mappings on the sphere, Proc. Amer. Math. Soc., 125(1997), 
51-60 

[12] I. P. Goulden and D. M. Jackson, Combinatorial Enumeration, Dover Publi¬ 
cations, Incorporated 2004, ISBN:0486435970 

[13] C. Itzykson, J. B. Zuber, Combinatorics of the modular group. 2. The Kont- 
sevich integrals, Int. J. Mod. Phys. A7 (1992) 5661-5705. |hep-th/920100lj. 

[14] M. Kazarian, KP hierarchy for Hodge integrals, Adv. Math. 221 (2009) 1-21. 
[arXiv:0809.3263 [math.AG]]. 

[15] M. Kazarian, S. Lando, An algebro-geometric proof of Wittens conjecture, J. 
Amer. Math. Soc. 20 (2007), no. 4, 1079C1089, [math. AG/0601760j. 

[16] M. Kontsevich, Intersection theory on the moduli space of curves and the 
matrix Airy function, Commun. Math. Phys. 147 (1992) 1-23. 

[17] W. Magnus, On the exponential solution of differential equations for a linear 
operator, Commun. Pure Appl. Math. 7 (1954), 649-673. 

[18] G. Marsaglia and J. C. W. Marsaglia, A New Derivation of Stirling's Approx¬ 
imation to n!, The American Mathematical Monthly, Vol. 97, No. 9 (Nov., 
1990), pp. 826-829, Published by: Mathematical Association of America, Sta¬ 
ble URL: http://www.jstor.org/stable/2324749 

[19] W. Miller, Symmetry Groups and their Applications, Academic Press, New 
York, 1972, pp 159-C161. ISBN 0-12-497460-0 

[20] T. Miwa, M. Jimbo, E. Date, Solitons: Differential Equations, Symmetries 
and Infinite Dimensional Algebras, Cambridge Tracts in Mathematics 135. 
Cambridge University Press, Cambridge, 2000 

[21] M. Mulase and B.Safnuk, Mirzakhani’s recursion relations, Virasoro con¬ 
straints and the KdV hierarchy, |arXiv:math/0601194 [math.QA]] (2006). 

[22] D. Mumford, Towards enumerative geometry on the moduli space of curves. 
In: Arithmetrics and Geometry (M. Artin, J. Tate eds.), v.2, Birkhauser, 
1983, 271-328. 

[23] A. Okounkov, Toda equations for Hurwitz numbers, Math. Res. Letters, 7 
(2000), 447-453, [arXiv:math/0004128j. 


52 



[24] M. Suzuki, Quantum analysis - Non-commutative differential and integral cal¬ 
culi. Communications in Mathematical Physics, 1997, Volume 183, Number 
2, Page 339. Stable URL: http://dx.doi.org/10.1007/BF02506410 

[25] V. S. Varadarajan, Lie Groups, Lie Algebras and Their Representations. 
GTM 102. Springer-Verlag, 1984. 


Xiaobo Liu 

Beijing International Center for Mathematical Research, 
Peking University, Beijing, China. 

E-mail address: xbliu@math.pku.edu.cn 

Gehao Wang 

Beijing International Center for Mathematical Research, 
Peking University, Beijing, China. 

E-mail address: gehao-wang@hotmail.com 


53 



